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Abstract 

Kink dynamics in spatially discrete nonlinear Klein-Gordon systems is considered. 
For special choices of the substrate potential, such systems support continuous trans- 
lation orbits of static kinks with no (classical) Peierls-Nabarro barrier. It is shown 
that these kinks experience, nevertheless, a lattice-periodic confining potential, due to 
purely quantum effects anaolgous to the Casimir effect of qunatum field theory. The 
resulting "quantum Peierls-Nabarro potential" may be calculated in the weak coupling 
approximation by a simple and computationally cheap numerical algorithm, which is 
C*~) ' applied, for purposes of illustration, to a certain two-parameter family of substrates. 
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Many systems in condensed matter and biophysics may be modelled by infinite chains of 
coupled anharmonic oscillators. If the anharmonic substrate potential has two or more 
degenerate vacua, such a system may support static kink solutions interpolating between 
neighbouring vacua. These kinks have various interesting physical interpretations (as crystal 
dislocations Q, charge density waves ||] and magnetic Q and ferroelectric [Q] domain walls, 
for example) and their dynamics is an interesting and important subject. 

Such a system of oscillators has an alternative interpretation as a spatially discrete 
version of an appropriate nonlinear Klein-Gordon equation. The lattice spacing h is related 
to the spring constant of the chain a by a = 1/h? so the strong spring-coupling limit 
is interpreted as a continuum limit. In the continuum limit, static kinks may occupy any 
position in space, by translation symmetry. This is generically untrue in the discrete system: 
static kinks may generically occupy only two positions relative to the lattice, one of which 
is a saddle point of potential energy, the other a local minimum. The difference in energy 
between these two static solutions is the Peierls-Nabarro (PN) barrier, the barrier which a 
kink must surmount in order to propagate from one lattice cell to the next. It can have 
strong effects on the dynamics of kinks in the system (kink trapping, radiative deceleration, 
phonon bursts etc. ||, §). 

One might expect the PN barrier, and hence its effects, to grow monotonically with 
h (of course the barrier vanishes as h — * 0). However, recent work of Flach, Kladko and 
Zolotaryuk JjJ has shown that this is certainly not universally true. In fact, there exist 
infinitely many substrate potentials with the property that at at least one none-zero lattice 
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spacing, h* say, the PN barrier vanishes exactly and a continuous translation orbit of static 
kinks is recovered. We shall say that a substrate potential with this property is "transparent 
at lattice spacing h*." Such potentials may be constructed by means of the so-called Inverse 
Method, and are clearly of some theoretical interest. 

The purpose of the present paper is to argue that although the kinks of such a sys- 
tem (at h = h*) are free of the classical PN barrier, they still experience a qualitatively 
similar periodic confining potential due to quantum effects analogous to the Casimir effect 
of quantum electrodynamics. We call this the quantum Peierls-Nabarro (QPN) potential. 
The essential physical observation is that the total zero-point energy of the lattice phonon 
modes in the presence of a kink depends periodically on the kink position. It should be 
emphasized that the kink position itself is treated as a classical degree of freedom while the 
phonons are quantized. The physical regime in which this is consistent will be identified: 
the classical kink mass must far exceed the phonon mass, and the kink must interpolate 
between widely separated vacua. For purposes of illustration, we shall compute the QPN 
potential numerically for a two parameter family of substrate potentials which (in a sense) 
includes discrete sine-Gordon and 4> systems. As a by-product of these calculations, we 
will obtain numerical evidence in favour of the assumption that kinks in these models are 
classically stable. 

2 Construction of transparent substrate potentials by the 
Inverse Method 

The general discrete nonlinear Klein-Gordon system consists of a field <p : Z x R — > R whose 
evolution is determined by a second order differential difference equation, 

4>n = J^(4>n+l -24>n + 4>n-l) -V'((j) n ). (1) 

Here h is the spatial lattice spacing, <\> n = d 2 (j) n /dt 2 and V is the substrate potential. One 
interpretation of the equation of motion is as that of an infinite system of identical oscillators 
(each oscillating in potential well V) with nearest neighbours coupled by identical Hooke's 
law springs of strength a = hr 2 . As the name suggests, ([lj) becomes a nonlinear Klein- 
Gordon equation in the continuum limit, h — > 0. If V(4>) has neighbouring degenerate vacua 
at (j> = a_ and <fi = a + > a_, say, (||) supports static kinks interpolating between them. To 
find these requires the solution of a second order nonlinear difference equation subject to 
the boundary conditions lim n _ »±oo 4> n = a±, which is usual only possible numerically. 

In this section we will construct, for a given lattice spacing h* > 0, a substrate potential 
Vh t (</>) which supports a continuous translation orbit of static kinks, and so by construction 
is transparent at lattice spacing h*. To do this we shall use a variant of the Inverse Method 
of Flach, Kladko and Zolotaryuk j7], |8| (which was originally devised to construct substrate 
potentials which support exact propagating, rather than static, kink solutions). The idea is 
to choose a static kink profile, that is an analytic, monotonic surjection / : R — > (a_,a_|_), 
satisfying exponential decay critera as z — > ±00, 

/(z) = o± + 0( e -"W), ^>0, (2) 
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and the symmetry requirement, 



f( z ) + f(-z) = a+ + a„ 



(3) 



then impose that the translated kink <\f n = f(nh — b) be a static solution of the system (with 
h = h*) for all & £ I. From equation (||), this condition holds provided one chooses V^, 
such that 

K(f(*)) = ^lf( z + h *) ~ 2 /(*) + /(* ~ Ml ( 4 ) 

for all z S R. This uniquely determines : (a_,a+) — > R by monotonicity of /, and 
hence Vh„ ■ (a_,o+) — > R up to an arbitrary constant. To complete the definition of this 
transparent substrate, one should extend its definition appropriately to all R. How one 
does this is somewhat arbitrary, but will have no bearing on our results, so we shall merely 
demand that Vh t , , be continuous at a±. Equations (0) and @ then imply 
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KM 



o 



hi 



(cosh [xh* 



1) > 0. 



(5) 
(6) 



Hence 4> = a± are both stable equilibria. To see that these are degenerate vacua, note that 
from (|), 



hl\V hm (a + )-V h .(a-)] 



/oo 
[f(z + K)-2f{z)+f{z-K)]f'(z) dz 
-oo 

/oo 
[f(z + K) + f(~(z + K))]f'(z) dz - (a< 
-oo 

(al - a 2 _) 



(a+ + a_)[/(*)] 
0, 



oo 
— oo 



(7) 



using the symmetry constraint on /, (^). 

So given a kink profile /, the inverse method generates a one-parameter family of double- 
well substrate potentials {Vh„ ■ h* > 0} each of which is transparent at spacing h = h*. 
Moreover, one has as explicit formula for the continuous translation orbit of static kink 
solutions, namely (f) b n = f(nh* — b), b G R. On physical grounds, one expects the static 
kinks to be stable to small perturbations, although strictly speaking this is not assured. 
One would need to check that the Hessian of the potential energy functional, 



Ep = j: 
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+ V h J 



(8) 



about 4> b has strictly positive spectrum (except for the zero-mode associated with transla- 
tion). The quantum calculation described in section |3] may be reinterpreted as the calcu- 
lation of this spectrum. The results of section || then constitute numerical confirmation of 
kink stability for the family of transparent substrates considered therein. 
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3 The quantum Peierls-Nabarro potential 



In this section we shall quantize the system using a weak coupling approximation, essentially 
following the method outlined in ||, adapted to the infinite lattice. The method has 
previously been applied in the spatially discrete context to a certain nonstandard lattice 
sine-Gordon model pC|] . 

One may regard Ep as a potential energy function on the inifinite dimensional space Q 
of sequences <ft : Z —* K. The vacua <p = a± lie at the bottom of identical potential wells, 
cut off from each other, and all configurations with kink boundary behaviour, by an infinite 
energy barrier. Assuming that V is transparent at the lattice spacing under consideration 
(e.g. V = Vh, and h = h*, as in section |2|), the continuous kink translation orbit is an 
equipotential curve in Q: the classical energy of a static kink is independent of its position. 
If the kinks are stable, this is a level valley bottom winding through Q. 

Quantum mechanically, a particle cannot sit at the bottom of a potential well, or on the 
low dimensional floor of a valley: it always possesses a zero-point energy dependent on the 
shape of the well bottom. In this section we will semi-classically quantize motion both in 
the vacuum and kink sectors of the system. A physical regime will be identified in which 
the kink is very heavy, so that the kink position b may be treated as a classical degree of 
freedom, while the comparatively light phonon modes orthogonal to the translation mode 
are quantized perturbatively, by Taylor expansion of Ep. Computation of the kink ground 
state energy then amounts to summing the zero-point energies of an infinite system of 
harmonic oscillators, resulting in a divergent series. The quantity of physical significance 
is not this energy, but rather the difference between the kink and vacuum ground state 
energies, which is expected to be finite. Since translation is not a symmetry of the discrete 
system, there is no reason to expect this energy to be independent of b, that is, one expects 
the quantum kink energy to vary periodically with kink position, which is the origin of the 
quantum PN potential. 

It is convenient to introduce a dimensionless coupling constant A into the model so that 
the (classical) Hamiltonian of the system is 



(9) 



where n n = 4> n is the momentum conjugate to <j) n . Assuming that V is transparent at 
spacing h, this system supports a continuous translation orbit of static kinks 

^ X = jf(nh-b), b£R (10) 

interpolating between a_/A and a + /\. The classical energy of these kinks is independent 
of b and clearly scales with A as 

E P [4> b > x ] = ^Ep^ 1 ). (11) 

The physical regime of interest is that of small A where the kinks interpolate between 
widely separated vacua, by (|i~0|), and are very heavy, by (|TT1) . In this regime, one may 
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approximate motion about any stable static solution <f> = <p/X by using a truncated Taylor 
series approximation for Ep[<j) + 5<p]: 

E P $+5<j>] = ^E P [<p] + W nmHn^m + O(A) (12) 

n,m 

where 

W nm = 5 nm [2 + h 2 V"(ip)] - 5 n , m+ i - <5 n ,m-i- (13) 
Since W is a real, symmetric matrix, there exists an orthogonal transformation R such that 

W nm = J2 R lkUklRlm (14) 

k,l 

with U diagonal. The diagonal entries of U are the eigenvalues A n of W, none of which 
is negative provided 4> is stable, as we are assuming. Introducing normal coordinates ^ n — 
Ylm Rnmdfim which have conjugate momenta r\ n = Ylm Rnm^m, the Hamiltonian for motion 
about <f> reduces to 

H = ±E P [p] + \ + + 0(A). (15) 

n 

Neglecting the O(A) remainder, this is the Hamiltonian for an inifinite set of decoupled 
harmonic oscillators of natural frequencies hr 1 ^/ A n . 

We now quantize in standard canonical fashion in the cases where <j> = a±/X (the 
vacuum) and 4> = 4> b ' X (the kink located at b). Let the W-matrices in these cases be 
denoted W vac and W K (b) repectively, with spectra {A^ ac } and {A^(b)}. In each case, the 
quantum correction to the ground state energy is the sum of the zero-point energies of the 
oscillators, 

n 

in units where h = 1 (recall that h denotes the lattice spacing of the system). In the case of 
kinks, one should omit from this sum the eigenvalue associated with the translation mode 
since, the kink being very heavy, this mode is treated classically. Actually this makes no 
difference since the corresponding eigenvalue vanishes, so one might as well sum over all 
modes, including the zero mode. 

Of course, the series (16) is divergent in both the vacuum and kink sectors, and must 



be suitably regulated and renormalized. To this end, we truncate the lattice symmetrically 
about the kink centre (so —no < n < tiq, assuming b G [0, h)) and consider the spec- 
tra {A n (N) : n = —l...,N} of the truncated W-matrices of order N = 2uq + 1. The 
renormalized ground state energy is then 



E P [(P ^ 



N 



m = - J ^ r J - + iY 1 ™E[Vl A « W*)! - VAT(A0L (17) 



n=l 



which one expects to be finite, given the exponential spatial localization of the kink (the 
large \n\ entries of the matrix W are essentially identical to those of W vac ). The finite 
size of the lattice perturbs the translation zero mode away from zero slightly, so one of 
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the kink eigenvalues A n (N, b) may be negative for some (N, b) (although it must vanish as 
N — > oo). This is why we have introduced an absolute value into equation (p7|), so that £{b) 
is the limit of a real valued sequence. In the limit N — ► oo lattice translation symmetry is 
recovered, so £(b) should be periodic with period h. 



4 The cosh ^ kink family 

If we now choose V = Vh, for some kink profile /, we see from equations ( |l3|) and (0) that 
the QPN potential depends on / only through /', because 



f'(z + h) + f'(z - h) 

/'(*) 



For purposes of illustration we will consider the one-parameter family 

1 



cosh^ z 



(18) 



(19) 



with fj, > 0. Note that this includes the cases of the sine-Gordon (f(z) = 2 tan -1 e z ) and (ft 4 
(f(z) = tanhz) kink profiles: fi = 1 and [i = 2 respectively. The corresponding transparent 
substrates in these two cases are shown in figure 1. 




Figure 1: Transparent substrate potentials generated by the sine- Gordon (fi = 1, left) and 
(ft 4 (fi = 2, right) kink profiles, with lattice spacings from h = 0.5 (bottom) to h = 3.0 (top) 
in steps of 0.5. 



The neighbouring vacua a± may be any real numbers separated by 

f°° dx 

"+-°- = L^- (20) 

Since f(z) has order e~^ exponential decay, V^'(a±) is given by equation (||), and the 
vacuum Il^-matrix takes the simple form 

= 2 cosh nh6 nm - (<5 n ,m-l + $n,m+l)- (21) 
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The spectrum of the system truncated to N lattice sites is easily computed: 

AT(iV) = 2(cosh y,h - 1) + 4 sin 2 ( ^W ) > n = 1, 2, . . . , iV. (22) 

In the limit JV — > oo, the spectrum densely fills the interval [2(cosh yh — 1), 2(cosh fih + 1)]. 
The eigenvalue problem for the kink W-matrix, 

xx/K /ia f'(nh + h-b) + f(nh-h-b) x 

W nm {b) = Jlfnh - b) ~ (Avm-l + <W+l), (23) 

is intractable analytically, and will be solved numerically in section ||. One can show, 
however, that the quantum kink energy is (for all b) lower than the classical kink energy, 
that is, the quantum energy correction is negative. To see this, let A(6) be the real diagonal 
matrix 

A nm (b) = h 2 5 nm [V£(f(nh - b)) - Vl'(a ± )}, (24) 

so that W K (b) = W vac + A(6). Let {A^(6)}, {A v n ac } and {T n (b)} be the eigenvalues of 
W K (b), W vac and A(6), each spectrum arranged in nonincreasing order (T\ > T2 > T3 > 
. . .). Then standard matrix perturbation theory [O] asserts that 

A^(6)<Ar + r x (6) (25) 

for all n, where Ti is the greatest eigenvalue of A(6), 

1^(6) = maxh 2 [Vl'(f(nh - b)) - V^'(o±)]. (26) 

n 

So if VJl(f(z)) < VJl(a±) for all z e R then T„ < and the perturbation of W vac by 
A(6) must reduce each eigenvalue A^(6) relative to its vacuum counterpart, with the result 
(always assuming kink stability) that 



£ty|A*(6)| - VKT) < 0- (27) 

n 

This condition on Vl'^) (maximum second derivative at the vacua) is quite natural, 
and can easily be shown to hold for the whole cosh _Al family (/x > 0, h > 0). Note that 

h 2 K{f{z)) - Vl'(a ± )] = g(z) - lim g(x), (28) 

x — >oo 

where 

(29) 



(30) 



Now g is differentiable and even, and 

g'(z) = — /isinh/icosh^ -1 z 



cosh z 




cosh z 


cosh(z + h) 


+ 


cosh(z — h) 



coslr" +1 (z + h) coslr" +1 (z - h) 



so the only critical point of g is z = 0, a local minimum, whence it follows that 5(2) < 
\\m x ^ OQ g{x) for all 2. 
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5 Numerical results 



Since the truncated kink VF-matrix W£ (b) is real, symmetric and tridiagonal, its eigenvalue 
problem is particularly easy to solve numerically. In this section we present data generated 
by implementing the QL decomposition algorithm for tridiagonal matrices with implicit 
eigenvalue shifts, outlined in [12]. 

The first thing to check is that, as expected, the spectrum of W K {b) is positive semi- 
definite with nondegenerate eigenvalue zero. The least and next-to-least eigenvalues of 
(b) for N odd, 3 < N < 90 were computed for a large sample of parameter values in 
the range 1 < fj, < 3, 0.5 < h < 10, < b/h < 0.5. The results were similar at all points 
sampled: the least eignevalue converges to 0, while the next-to-least converges to a positive 
number below the lower edge of the vacuum phonon band, that is, less than 2(cosh [ih — 1). 
It is instructive to look at the build up of the spectrum of (b) as N grows large, as 
depicted for two contrasting sets of parameter values in figure 2. This clearly shows the 
rapid convergence of the bottom eigenvalue to (convergence being faster for larger \xh, since 
the kink structure is then more tightly spatially localized) and the next lowest eigenvalue 
to a constant outside the phonon band. The rest of the eigenvalues accumulate, apparently 
densely, within the phonon band (delimited by horizontal dashed lines in figure 2). This is 
precisely the right behaviour to ensure both kink stability and convergence of the quantum 
corrected kink energy in the limit N — > oo. The size of truncated system needed to obtain 
practical convergence depends on /j,h, but, within the parameter range cited above, N = 51 
seems to suffice. This is the matrix order used to obtain the remaining numerical data. 



■ iiii iii liil II II 



10 20 30 40 50 60 70 



90 100 



cn 3 cp j cm i cm i cm i cp t cm c cm c cm cm a cc 



Figure 2: Build up of the spectrum of (&) as the matrix order N grows large, in the 
cases \i = 1.5, h = 0.5, b/h = 0.2 (left) and \i = 2.0, h = 2.6, b/h = (right). Note how 
the spectrum accumulates within the vacuum phonon band, delimited by horizontal dashed 
lines. 



The quantity £{b) as defined in (|T^) is problematic to plot since jt contains contributions 
of different orders in A. For this reason it is convenient to consider £{b) = £(b) — £(0), which 
is of order A , instead. In fact, £{b) is the quantity of most direct physical significance 
anyway: it gives the change in kink energy as b varies, which is precisely what is meant by 
the quantum Peierls-Nabarro potential. 
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Figure 3 shows £{b) (0 < b/h < 1) for the sine-Gordon substrate potentials (/x = 1), at 
a variety of lattice spacings. The results are qualitatively very similar to the usual classical 
PN potential: the kink has greatest energy when located exactly on a lattice site (b = 0) 
and least when located exactly midway between lattice sites (b = h/2), the energy difference 
(the QPN barrier) growing monotonically with h. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 3: Position dependence of the quantum PN potential £{b) for fi = 1 and h = 1.5 to 
h = 5 in steps of 0.5. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 

b/h 



Figure 4'- Position dependence of the quantum PN potential £(b) for fi = 2 and h = 2 to 
h = 5 in steps of 1. The unlabeled curves are h = 2 (dashed) and h = 4 (solid). 



Figure 4 shows similar plots for the 4> 4 substrate potentials (// = 2). Here, once again, 
kinks have greatest energy when 6 = and least when b = h/2, but the QPN barrier does 
not grow monotonically with h, nor is the shape of £(b) so uniform as in the fj, = 1 case. 
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Similar plots for \i = 3 reveal more complicated behaviour, as shown in figure 5. For 
small h, £{b) is maximum at b = and minimum at b = h/2 as for // = 1, /i = 2. However, 
for h above a critical value (~ 1.46) 6 = becomes a local minimum of £ and two extra local 
maxima appear. The global minimum of £ remains at b = h/2, rather than 6 = 0, until h 
exceeds about 1.52, after which £(h/2) exceeds £ (0). As h is increased further, the two local 
minima coallesce at 6 = h/2 (at h ~ 1-7), so that the QPN potential starts to resemble an 
inverted version of the \i = 1 case: now kinks have minimum energy when located exactly 
on a lattice site and maximum energy when located exactly midway between sites. 



0.025 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

b/h 

Figure 5: Position dependence of the quantum PN potential £(b) for [i = 3 and h = 1.4 to 
h = 1.6 in steps of 0.2. 




3 4 5 6 7 8 9 



Figure 6: A rough measure of the depth of the QPN barrier: £(h/2) as a function of h for 
various fi G [1, 3]. 
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Periodicity and reflexion symmetry of £ imply that b = and b = h/2 must always be 
critical points. There may be others (as in the case fi = 3.0), but for the most part, plotting 
£(h/2) = £(h/2) — £(0) against h gives a good account of how the QPN barrier varies 
with h and (i. In particular, the sign of £{h/2) tells one whether the QPN barrier tends 
to trap kinks between lattice sites (6 (h/2) < 0) or on lattice sites (£(h/2) > 0). Figure 
6 shows plots of £(h/2) against h for various jj, G [1,3]. Three regimes clearly emerge for 
large h: 1 < [i < 2 (£(h/2) < growing unbounded as h — ► oo), \i = 2 (£(h/2) < 
remaining bounded as h — > oo) and 2 < /u < 3 (£(h/2) > growing unbounded as h — > oo). 
This trichotomy may be explained by consideration of the asymptotic forms of W K (0) and 
W K (h/2) for large h: 

W K (0) 



W K (h/2) 

e fJ,h 



diag(...,l,l,2^,0,2^,l,l,...) 
diag(...,l,l,0,0,l,l,...). 



This leads to the prediction 



£(h/2) 

g—/j,h/2 

• £(h/2) 



1 

2h I 
(1 



(0 - 0) + (0 



-A*/2> 



+ (1 



-A*/2> 



+ (1-1) + (1-1) + 



fih/2 



2h 



(31) 
(32) 

(33) 
(34) 



Formula ( |34"| ) accounts well for the asymptotic behaviour seen in figure 6. Clearly the most 
interesting case from this point of view is the critical value fj, = 2. The dependence oi£(h/2) 
on h for /j, = 2 is shown in figure 7. One sees that, rather counterintuitively, the QPN barrier 
actually vanishes in the extreme discrete limit, h — > oo. One should remember, of course, 
that in varying h one also varies the shape of the transparent substrate potential Vh (which 
would not otherwise remain transparent). In fact, the limit h — > oo is always badly singular 
since the curvature of the substrate at the vacua grows unbounded, by equation (||). 




Figure 7: Dependence of the QPN barrier on the lattice spacing h in the critical case (fx 
2). 
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6 Concluding remarks 



In this paper we have considered oscillator chains with no classical PN barrier and shown 
that their kinks still experience a lattice-periodic confining potential due to purely quantum 
mechanical effects, leading to a new mechanism for kink pinning. The quantum PN potential 
was computed numerically for a simple two-parameter family of substrates, revealing a rich 
variety of behaviour. 

It remains to be seen whether the QPN potential has any relevance to genuine physics. 
Given the idealized one-dimensional nature of the model, this seems unlikely (generalizing 
the Inverse Method to higher dimensions is very problematic). Another cause for doubt is 
that the effect only exists for certain special substrates. Just how special these "transparent" 
substrates are is unknown. In order to have any physical relevance they would at least have 
to be structurally stable: if Vh[f] is the substrate transparent at h generated by kink /, 
then given any sufficiently small perturbation 5V there should exist a kink /* close to / and 
spacing h* close to h such that Vh[f] + 5V = Vh, [/*]. Thinking of the Inverse Method as a 
mapping K x R + — > P (where K is the space of kink profiles and P the space of potentials), 
the question is whether this mapping is continuous with respect to some sensible choice of 
topologies on K and P. 

This is one of many open mathematical questions raised by the present work. We have 
presented numerical evidence to support the assumption of classical kink stability, but it 
should be possible to prove stability rigorously. Similarly, one should be able to prove 
convergence of the series defining £ (b). In both cases one needs to understand the large ./V 
behaviour of the spectrum of W^(b). Standard minimax estimates are insufficient for this 
purpose - more delicate analysis is required. 

Acknowledgments 

This work was partially completed at the Max-Planck-Institut fur Mathematik in den Natur- 
wissenschaften, Leipzig, where the author was a guest scientist. He wishes to thank Prof. 
Eberhard Zeidler for the generous hospitality of the institute. The author is an EPSRC 
Postdoctoral Research Fellow in Mathematics. 

References 

[1] Y.Y. Earme and J.H. Weiner, "Loss-free dislocation motion in a lattice model" Phys. 
Rev. Lett. 33 (1974) 1550. 

[2] M.J. Rice, A.R. Bishop, J. A. Krumhansl and S.E. Trullinger, "Weakly pinned Frohlich 
charge-density-wave condensates: a new, nonlinear, current-carrying elementary exci- 
tation" Phys. Rev. Lett. 36 (1976) 432. 

[3] A.R. Bishop and W.F. Lewis, "A theory of intrinsic coercivity in narrow magnetic 
domain walls" J. Phys. C12 (1979) 3811. 

[4] D.A. Bruce, "Scattering properties of ferroelectric domain walls" J. Phys. C14 (1981) 
5195. 



12 



[5] M. Peyrard and M.D. Kruskal, "Kink dynamics in the highly discrete sine-Gordon 
system" Physica 14D (1984) 88. 

[6] R. Boesch, C.R. Willis and M. El-Batanouny, "Spontaneous emission of radiation from 
a discrete sine-Gordon kink" Phys. Rev. B40 (1989) 2284. 

[7] S. Flach, Y. Zolotaryuk and K. Kladko, "Moving lattice kinks and pulses: an inverse 
method" Phys. Rev. E59 (1999) 6105. 

[8] J.M. Speight, "Topological discrete kinks" Nonlinearity 12 (1999) 1373. 

[9] R. Rajaraman, Solitons and Instantons (North-Holland, Amsterdam, 1989) Ch. 5. 

[10] J.M. Speight, "Kink Casimir energy in a lattice sine-Gordon model" Phys. Rev. D49 
(1994) 6914. 

[11] J.H. Wilkinson, The Algebraic Eigenvalue Problem (Oxford University Press, Oxford, 
1965) pp 101-102. 

[12] W.H. Press, B.P. Flannery, S.A. Teukolsky and W.T. Vetterling, Numerical Recipes 
(Fortran Version) (Cambridge University Press, Cambridge, 1989) pp 362-363. 



13 



